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Abstract. We study magnetically powered relativistic outflows in which a part of the magnetic energy is dissipated 
internally by reconnection. For GRB parameters, and assuming that the reconnection speed scales with the Alfven 
speed, significant dissipation can take place both inside and outside the photosphere of the flow. The process leads 
to a steady increase of the flow Lorentz factor with radius. With an analytic model we show how the efficiency of 
this process depends on GRB parameters. Estimates are given for the thermal and non-thermal radiation expected 
to be emitted from the photosphere and the optically thin part of the flow respectively. 

A critical parameter of the model is the ratio of Poynting flux to kinetic energy flux at some initial radius of the 
flow. For a large value 100) the non-thermal radiation dominates over the thermal component. If the ratio is 
small 40) only prompt thermal emission is expected which can be identified with X-ray flashes. 
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1. Introduction 

To overcome the compactness problem of 7 -ray bursts 
(GRBs) (e.g. Piran 199^ ) the central engines must pro¬ 
duce radiating material moving ultra-relativistically fast 
towards the observer. GRB models must therefore de¬ 
scribe an energy source which not only releases energy of 
around 10 ®^ erg/sterad but must also explain the ‘clean’ 
form of the energy. To produce the high Lorentz factors 
of the orde r of 10^-10^ ([Fenimore et al. 1993 ; Woods & 


Kluzniak fc Ruderman 1998| ; [Spruit 19991) would produce 
such a rotationally driven Poynting flux. Extraction of en¬ 
ergy from the central object by this magnetic process is 
potentially very efficient and fast. 

In order to obtain not only a large energy extrac¬ 
tion but also the observed large bulk Lorentz factors, the 
Poynting flux must be converted to kinetic energy. The 
simplest available magnetic acceleration models, in which 
the flow is approximated as radial, are problematic in 


[Loeb 1|995|; Ewick fc Sari 200l|) which are needed only non-relativistic case ( jWeber 

^ fc Davis 1967| ; [Belcher fc MacGregor 197^ a dominat- 


a sm rf fraction of the total en e rg}^ can exi s t in form of 
rest mass energy of the matter involved. 

The popular models involving compact objects or the 
collapse of a massive star to a black hole must include 
mechanisms how the energy is transported into a space 
region with few baryons. Otherwise large amounts of mass 
are expelled which cannot be accelerated to high Lorentz 
factors. The initially released energy could leave the cen¬ 
tral polluted region by ne utrinos which annihilate to a 
pair plasma further away (Berezinskii fc Prilutskii 1987; 
Goodman et al. 1987 ; Ruffert et al. 1997 ). But due to the 
small cross section of neutrinos the efficiency is low and 
most o f the energy oscapoB as neutrinos. - 


ing initial Poynting flux can transfer 1/3 of its energy to 
the matter. If the flow is initiall y relativistic however al¬ 
most no acceleration is possible ( Michel 1969 ). The phys¬ 
ical reason lies in the singular held and flow geometry of 
a purely radial flow. In this case the magnetic pressure 
gradient balances the magnetic tension force and no ac¬ 
celeration occurs. An imbalance between the pressure gra¬ 
dient and the tension force occurs in non-radial outflows, if 
the flow lines diverge faster with radius than in the radial 
case (Begelman fc Li 1994; Takahashi fc Shibata 1998). 
Detailed 1-dimensional calculations have been made which 


show how such a flow divergence can come about (Beskiu 


A Poynting flux dominating outflow will naturally oc- 
cur if the compact object rotates and possesses a magnetic 
held. The luminosity will be fed by the rotational energy 
reservoir of the central object. Models involving an mag¬ 
netised torus around a black hole ( Meszaros fc Rees 1997 ) 
or a highly magnetised millisecond pulsar (Usov 1992; 


1997; Daigne fc Drenkhahn 200^ ). 
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In this paper we show that there is a second process 
which naturally leads to efficient conversion of Poynting 
flux to bulk kinetic energy. If the magnetic held in the 
outflow contains changes of direction on sufficiently small 
scales, (a part of) the magnetic energy is ‘free energy’ 
which can be released locally in the flow by ‘fast recon¬ 
nection’ processes. Such a decay of magnetic energy, if it 


























































2 


G. Drenkhahn: Acceleration of GRB outflows by Poynting flux dissipation 


can occur rapidly enough, has two desirable effects. First 
it provides a source of energy outside the photosphere 
which is converted directly into radiation, without the 
relatively inefhcient intermediate step of internal shocks 
( ^pruit et al. 2001 , hereafter Paper I). Secondly, it leads to 
an outward decrease of magnetic pressure, which causes a 
strong acceleration of the flow and conversion of Poynting 
flux to kinetic energy. In the present work, we concentrate 
on the acceleration effect, and show how it depends on the 
parameters (energy flux, baryon loading) of a GRB. This 
aspect of the model can be illustrated with analytic cal¬ 
culations. In a future paper we show, with more detailed 
numerical results, how the dissipated magnetic energy can 
also power the observed prompt radiation. 

Changes of direction of field lines must occur in the 
flow in order for energy release by reconnection to be 
possible. These can occur naturally in a number of ways. 
If the magnetic field of a rotating central object is non- 
axisymmetric the azimuthal part of the magnetic field in 
the flow changes direction on a length scale A ~ ttu/H, 
where v is the flow velocity and If the angular frequency. 
For an inclined dipole this yields the ‘striped’ field in pul¬ 
sar wind model of Coroniti (1990) where magnetic energy 
is released by the annihilation of the antiparallel field com¬ 
ponents. Field decay by reconnection was applied to pulsar 
winds dCoroniti I990|; Lyubarsky & Kirk 2001) and also 
to GRBs (Thompson 1994; Paper ]| )! 

In this paper we investigate the dynamics of a magnet¬ 
ically powered outflow in which some of the energy dissi¬ 
pates by reconnection. With the assumption that the flow 
is highly dominated by magnetic energy and that the ther¬ 
mal energy is negligible we derive the velocity profile of 
the flow. The results provide estimates of the Lorentz fac¬ 
tor of the flow, the photospheric radius, and the amount of 
energy that can be converted into non-thermal radiation. 
We investigate under which conditions prompt emission is 
expected and whether a considerable amount of thermal 
radiation can be produced. These predictions can then be 
tested against observations of the thermal component in 
GRB spectra (Preece 200C). 


2. Model description 

Highly magnetised spinning compact objects, e.g. mil¬ 
lisecond pulsars or tori around black holes, are sources 
of Poynting flux that can power GRBs. They produce 
a plasma-loaded electromagnetic wind travelling outward 
and are fed by the rotational energy of the central ob¬ 
ject. In the wind of an aligned rotator the magnetic field 
is ordered and stationary. If ideal MHD applies, and the 
wind is radial in the poloidal plane, a large fraction of 
the total luminosity is bound to stay in form of Poynting 
flux. The picture changes in the case of an inclined ro¬ 
tator or any other source producing a non-axisymmetric 
rotating magnetic field. If the emitted Poynting flux con¬ 
tains modulations of the field it also carries along free 
magnetic energy, which can be extracted by reconnection 
processes. In these processes the field rearranges itself to a 


energetically preferred configuration while the energy re¬ 
leased is transfered to the matter. Because perfect align¬ 
ment of magnetic and rotation axis is a special case it is 
likely that most astrophysical objects produce modulated 
Poynting fluxes containing free magnetic energy. 

A necessary condition for the existence of free magnetic 
energy in the flow is the field variation on small scales. 
For reconnection processes differently oriented field lines 
must come close to each other. Therefore the length scale 
on which the orientation of magnetic field lines change 
controls the speed of the field dissipation. The smaller the 
length scale is the faster the field can decay. 

The general large scale magnetic field structure ex- 
pe cted to b e produced by a rotating object was discussed 
in Paper I. It is useful to consider simplified flow geome¬ 
tries along the equatorial plane and along the rotation axis 
as examples. In the equatorial plane an inclined rotator 
will produce a ‘striped’ wind ( [Coroniti 1990 ). It consists 
of an electromagnetic wave in which the azimuthal field 
component varies with a wave length of 27rc/H. Along the 
rotation axis the wave will have a circular component with 
the same wave length. Such wave-like field variations are 
present in general if a non-axisymmetric magnetic field 
component is present. The equatorial plane of an inclined 
rotator is only a prototype to illustrate the field geometry. 
In general, wave-like variations occur at all latitudes. If the 
rotator is aligned the field will be axisymmetric. Then, the 
magnetic field geometry looks like a wound up spiral on 
all cones of equal latitude. This field geometry is present 
in case of a jet-like outflow. Here, the magnetic field does 
not vary on small scales along the outflow direction. The 
differently directed field components lie on opposing sides 
of the rotation axis. In the context of a jet-like outflow the 
typical length scale of the field variation is the diameter 
of the jet cone rd where d is the jet opening angle. 

In both of these field geometries MHD instabilities can 
promote reconnection processes. For wave-like variations 
current sheets form and tearing instability will lead to 
reconnection. For a polar jet-like outflow of an aligned 
rotator the field con figuration is hig hly unsta ble to th e 
kink instability (e.g. Bateman 1980| , see also Paper I ). 
It is plausible that the kink instability working in this 
case will distort the geometry after some time so that also 
wave-like variations come into play. This leads to non¬ 
periodic and highly irregular waves. Our model assumes 
the longitudinal field variation to be periodic so that the 
complicated effects of any non-periodicity is neglected. We 
consider both limiting cases for the small scale field varia¬ 
tions though wave-like structures seem to be more general. 

Near the source the flow is accelerated magnetocen- 
trifugally (and perhaps thermally). It will be accelerated 
up to a distance around the Alfven radius and then start to 
become radial asymptotically. The poloidal and azimuthal 
field components at the Alfven radius are similar in mag¬ 
nitude. Beyond this point their ratio scales as B^/Br ^ r, 
so that the radial component soon becomes negligible at a 
couple of Alfven radii. The Alfven point lies always inside 
the light radius c/H and if the magnetic field dominates, 
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like in our case, the Alfven radius and light radius are 
almost equal. Thus we can simplify the flow and fleld ge¬ 
ometry at source distances r > a few x c/Vl « 10^ cm 
by assuming a purely radial flow with an azimuthal mag¬ 
netic field. At this distance gravity effects can also be ne¬ 
glected. The magnetocentrifugal effects accelerate the flow 
to fast magneto-sonic velocity. Because we work in the cold 
limit (see below) and approximate the magnetic fleld to 
be purely poloidal the magneto-sonic velocity is equal to 
the Alfven velocity. The initial flow velocity is set to the 
Alfven velocity at some initial radius ro ^ c/0. 

To make a simple approach feasible analytically we 
have to make further approximations. The flow is treated 
stationary and its thermal energy is neglected (‘cold’ 
limit). This ‘cold’ approximation is quite good in the opti¬ 
cally thick region since no energy can be lost by radiation 
anyway. All of the dissipated energy is always converted 
into kinetic form. If the flow is optically thin the radia¬ 
tion produced by dissipation will freely escape and this 
energy part will not be converted into kinetic energy. Our 
model overestimates the kinetic energy gained in the op¬ 
tically thin regime. Statements about the radius of the 
photosphere, where the flow changes from optically thick 
to thin, or the Lorentz factor there will hold rather ro¬ 
bustly. 


2.1. Magnetic field dissipation 


The dissipation is modelled by using the typical length 
scale of the magnetic field A and a fraction of the Alfvenic 
velocity evA where e is an dimensionless factor. The idea 
is that the magnetic field lines with different directions 
get advected towards a reconnection centre where the field 
dissipates (Petschek 1964). This advection happens with a 
fraction of va- The decay of mean field must also depend 
on Va in a similar fashion. Though there are elaborate 
models on the reconnection physics (e.g. Coroniti 1990; 
Thompson 1994) we prefer to express this rather uncertain 
topic in form of the dimensionless free parameter e and the 
Alfvenic speed va- All the uncertain physics in this picture 
is taken up by e. Our ansatz for the time scale of the mean 
field decay, in a comoving frame is then 


Ac 


er'A,( 


( 1 ) 


where quantities considered in the comoving frame are 
indexed with ‘co’. This comoving frame moves with the 
mean large scale bulk flow motion so that the small scale 
motion is neglected. 

The reconnection takes place at certain reconnection 
centres in the flow. The typical distance between these 
reconnection centres also influence the rate of the overall 
field dissipation. Since we regard the field dissipation in 
all generality and for a variety of field geometries we can¬ 
not explicitely model these small scale details about the 
density of reconnection centres. This uncertain issue must 
also be handled by the free parameter e so that e controls 
the average field dissipation on larger length scales. 


At first sight e < 1 is an upper limit since for e = 1 
reconnection would happen everywhere with an advection 
speed of c. If the advection towards the reconnection cen¬ 
tres happens with almost c large current densities are re¬ 
quired. The MHD condition might break down leading to 
an additional decay of magnetic field. This effect could be 
parameterised by an larger value of e so that an upper 
limit of 1 may not be strict. 

For most of the paper we will work with a fiducial value 
of e = 0.1. One should keep in mind that e is perhaps the 
most uncertain quantity of the model because it may not 
be constant and its value cannot be estimated by first 
principles in general. 

The length scale for the dissipation Aco depends on the 
nature of the outflow as discussed in the last section. We 
will distinguish the two cases where the field variation is 
encountered on length scales longitudinal to the flow di¬ 
rection (called longitudinal case in the rest of the paper) 
and where this length scale is transversal to the flow di¬ 
rection (transversal case). The transversal case is found in 
a polar outflow of a aligned rotator where the field com¬ 
ponents having different directions lie on opposite sides 
of the rotation axis. For mathematical simplicity we will 
regard here the two limiting cases only and make a few 
notes on the mixed case later in this study. 

The longitudinal and the transversal length scales 
Aio, Atr in the comoving frame scale differently with the 
Lorentz factor of the flow F: 

Aio.co = ^ , (2) 

Atr.co = dr (3) 

where d is some kind of an opening angle of the polar 
outflow. Atr.co does not scales with F because it denotes a 
length which is perpendicular to the direction of motion 
being the same in the lab and comoving frame. 


3. The dynamics of the flow 

The reconnection processes described in the last section 
will change constantly the structure of the magnetic field 
on a length scale of the order of the wave length (small 
length scale of the problem). Though, the azimuthal field 
line stretching will keep the field aligned predominantly 
perpendicular to the flow direction. The exact field struc¬ 
ture is not important because only the magnetic energy 
density i?^/(87r) enters the dynamic equations. In the 
following, B denotes the dynamical effective transversal 
magnetic field which is constant over small scales. The 
induction equation will still be valid for the effective field. 

3.1. Conservation laws 

The dynamics of the flow is governed by the ideal MHD 
equations for the conservation of energy, momentum and 
mass. For the relativistic treatment the equations are for- 
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mally best written in tensorial form (e.g. Bekenstein & 
Oron l[^: 


= 0 , 
(pw'");/. = 0 


( 4 ) 

( 5 ) 


where 


'■ -.-^ 

matter part 

&a6“ - b^^b'' 
electromagnetic part 



( 6 ) 


is the energy-momentum tensor for ideal MHD. The sig¬ 
nature (—h -f-l-) is used for the metric tensor g^’'. Here, 
p, w and p are the mass density, the enthalpy density and 
the pressure in the proper frame of the fluid moving with a 
4-velocity = (r,M). b^ = *F^''ui, is the the 4-vector of 
the magnetic field where *F^'^ is the dual electromagnetic 
field strength tensor. 

Now we choose a spherical coordinate system centred 
on the central engine. The flow is assumed to be spher¬ 
ically symmetric and the held dominated by its toroidal 
component. In this case u F B and the components of 
the magnetic four vector are simply b^ = {0,—B/r) and 
bfj}^^ = Writing (||) and (||) in coordi¬ 

nate form and assuming stationarity gives the conserva¬ 
tion laws for energy, momentum and mass 


0 = dy(wru + !^^ , 

0 = oy + (l + + r^^rP , 

0 = dy pu 


( 7 ) 

( 8 ) 
( 9 ) 


where (3 = u/T ( [Konigl fc Granot 2002 ; Lyutikov 2001] ). 
By integrating the mass and energy equations one obtains 
the total mass loss per time per sterad 


M = r^puc 


( 10 ) 


and the total luminosity per sterad 


L = r^ 


v..r«c + y 

An 


( 11 ) 


Taking the flow to be cold and eliminating {tBY from 
(|T^, ( |T^ shows that it is a constant function of r. Finding 
an exact accelerating solution of the energy and momen¬ 
tum equations without thermal pressure is not possible. 
By using an evolution equation for the magnetic held B 
(see Sect. |3.4| below) and combining it with the energy 
equation (|13|) one obtains an accelerating solution but 
violates the momentum conservation. Luckily, the error 
made by that becomes small in the ultra-relativistic limit. 
Then, F « it, /3 « 1 so that the first term in the momen¬ 
tum Eq. (^) becomes small since it approaches the form 
of the rhs of the energy Eq. (|^). As a consequence the 
thermal pressure gradient term of (H) must also be small. 
Setting the pressure to zero and solving the energy equa¬ 
tion means that the momentum equation is almost sat¬ 
isfied. Since we only consider ultra-relativistic flows the 
error made is small which justihes the use of the cold ap¬ 
proximation. 

3.2. A note on the ideal MHD approximation 

In the treatment above the ideal MHD approximation was 
used. But a key ingredient of the model is the existence of 
field dissipation for which ideal MHD is not applicable at 
first sight. The held dissipation acts like an effective diffu- 
sivity in the plasma so that the effective mean electric held 
in a comoving frame does not vanish. Since a substantial 
electric held Eco would contribute to the comoving energy 
density, the question arises if it can be neglected. We found 
from a more detailed numerical investigation (in prepara¬ 
tion) that the comoving electric held is in fact small, and 
we use this advance knowledge to neglect its contribution 
to the dynamics here. 

3.3. The initial conditions of the flow 

Let (7 be the ratio of Poynting hux to matter energy flux: 

. = = = ,15) 

Akin dTrFMc AnT^pc^ 

a is also the magnetisation parameter of the plasma, de¬ 
scribing the ratio of the proper magnetic energy density 
to the proper energy density of the matter. The Alfven 
4-velocity in the comoving frame is 


The enthalpy density w includes the rest mass energy 
density pc^. In the following we will assume a cold flow 
with w = pc^, p = 0. Then, the momentum equation can 
be integrated and the conservation laws read 

M = r'^puc , (12) 

L = TMY + , (13) 

47r 

(tBV- 

const. = uMc + (1 + p3) - • (14) 

Stt 

In the energy equation (H) one can identify the kinetic 
energy flux per sterad Lkin = FMc^ and the Poynting 
luminosity per sterad Lpf = /3c(rB)^/(47r). 


UA = 


B/T 


\J Anpc^ 

with the regular dimensional velocity counterpart 
UA 


VA= c 




= C 


1 + cr 


(16) 


(17) 


At an initial radius Tq, where the flow starts with the 
Alfven speed (discussed in Sect. ||), the relation between 
initial 4-velocity uq and initial Poynting flux ratio (Tq is 
simply 


Mo = u{ro) = ua(?'o) = yp . 


(18) 
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The total energy and the mass flux are linked by 
L= {ao + l)ToMc^ (19) 

so that M can be expressed in terms of L and ctq: 

L 


M = 


c2(ao +1)3/2 


( 20 ) 


In the GRB case the flow must start highly Poynting 
flux dominated with (Tq ^ 100 so that m 1 at all dis¬ 
tances. One can therefore set co -I- 1 « cto, u ~ T and 
/3 « 1. The conservation Eqs. dl)-® then reduce to 

M = r'^puc , (21) 

/ \ {tBY 


L\l- 


3/2 


At: 


( 22 ) 


In this limit one can also simplify the expression (|lq ) for 
the Alfven speed in the comoving frame by using ( |2l|) and 


3/2 


Pa = Ji- 


3/2 


UA = 


(23) 


Equation (|^) states that there is no acceleration if 
i? ~ 1/r. This is the case in a radial outflow with ideal 
MHD conditions. We encounter here again the fact that 
the Poynting flux energy in a radial ultra-relativ istic MHD 
outflo w cannot be transfered to the ma tter (e.g Begelman 
& Li 1 Daigne fc Drenkhahn 2002 ). 


3.4. The evolution of the magnetic field 


The evolution of the magnetic field, as it is carried with 
the flow, is governed by the induction equation. This in¬ 
cludes the effects of advection and field line stretching. In 
addition we will include a term to describe the decay of (a 
part of) the field by reconnection as described in Sect. |2.l| . 
Since the reconnection is easiest described in a local, co¬ 
moving frame, we first transform the induction equation 

dtB = —ccxivlE (24) 

into the comoving frame where we extend it to account 
for the reconnection. 

In the stationary case of our model setup the induction 
equation for ideal MHD is 

drdrB = 0 . (25) 


This equation describes the field evolution due to ideal 
MHD processes. To obtain the evolution term in the the 
comoving frame we first need the convective derivative 

^ = cddrB = -—drrd (26) 

at r 

which then gives in terms of the comoving quantities 


^ ^ d(rHeo) 
dt Tdtco 

_ dHco Bco dP 

dtco 1/ dtco 

_ dHco B dr 

dtco r dt 


(27) 

(28) 
(29) 


Combining ( p^ ) and (^) gives the comoving held evolu¬ 
tion without dissipation effects 


dHco cB 

-J— = -—drTU . 


(30) 


Let us denote the striped, decayable part of the 
magnetic field with and the perpendicular, non¬ 

reconnecting part with B.^^. We model the decay of the 
striped component of the magnetic held in the comoving 
frame by the ansatz 


dtco 


cH- 

Pr 



(31) 


where r is the field decay time scale from (|^) in the lab 
frame. The non-decaying part evolves according to the 
induction Eq. (^ ) so that 


/3rH.|Y = const. 


(32) 


Expressing the como ving quantities in terms of lab frame 
quantities similar to (P^-(29) gives 


dH^co _ dB^ B^ dP 
dtco dt r dt 


(33) 


Since the flow is stationary we can replace the convective 
derivatives by r-derivatives and after combining (^ ) and 
( ® ) one obtains 

cS—^ B ^ 

cfddrB^ - T^cddrT = ——^drru -^ . (34) 

r Pr T 

One arrives at 

drdrB^ = (35) 

CT 


and with j3rB^ = sJ{j3rBY — {dcBY"^ this yields 


drfdrB 



(36) 


Let g = {B.^/B)q so that g is the initial fraction between 
the field strengths of the non-decaying component and 
the total field, g = 0 means that the complete field de¬ 
cays while at g = 1 only ideal MHD processes occur and 
no dissipation takes place. In the case of an equatorial 
outflow of an inclined rotator g = cost where i is the in¬ 
clination. The other cases are more complicated and g is 
not associated to a simple geometric quantity. 

Using /i as a constant of the problem the field decay 
equation can be cast into this simple form: 


dr{(3rBf = -2 


[fir by 
PCT 


_ APrBYp - 
^ {3rBY_ 


(37) 


Equations (|^), (p^), (p^, ( ^^ over-determine the 3 
unknown functions p, u, B because we have assumed a cold 
flow so that the internal energy is neglected. But in the ul¬ 
tra relativistic limit the energy and momentum Eqs. (p^, 
(®) are equal. In this limit (PH), d^), d^ ) are sufficient 
to obtain the solutions for p, u, B. 
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3.5. Solution of the flow problem 

In the ultra relativistic limit /3 « 1, F « u Eq. ( |^ ) for the 
evolution of the magnetic field with distance, including 
dissipation becomes 


dr-irB)^ = -2 


{rBf 

CT 



ArB)l - 

^ {rB)\ 


(38) 


Using the B-u relation to eliminate (rB)^ one obtains 
a differential equation for u: 



a-y^ (l - - u 


(39) 


This equation is to be integrated with initial condition 
u = uq = '/S. In the absence of internal dissipation {p = 
1) the flow is not accelerated, drU = 0, as expected. The 
flow accelerates monotonically and reaches asymptotically 
its terminal speed Moo found be setting drU = 0: 


Uoo = (l - p^) + ■ 

The dissipation time scales are 
27r V? 


T\o = 


eft 


\J 1 - u/ay 


(40) 


(41) 


and 


t) ru 

ec r , 3/2 


(42) 


for the longitudinal and transversal cases from (|^, (^, 

(i, (H- 

Since our model rests on the assumption of a significant 
‘decayable’ component, p in the following is taken to be 
of the order 0.5 but not close to 1. Then, the terminal 
velocity Uoo is much larger than the initial velocity and 
Uao/uo « cro(l — /r^) ^ 1 and (^) simplifies to 


^ /, 2\ 3/2 

Uoo « (1 - /i jCTo' . 


(43) 


The differential Eq. ( ^^ is analytically solvable at 
intermediate source distances, where the flow is much 
faster than the initial velocity (^/oo u) but is still 
far away from the point where the acceleration saturates 
(m <C Moo < rry^)- The dissipation time scales then sim¬ 
plify to 


2ttu'^ 


T\o 


Ttr 


ftru 


efl ec 

In this case Moo — m « Moo and 
2 ?; 


3) becomes 


drU = 


CT 


(44) 


(45) 


with the solutions 


—efluoo{r - ro) -I- ) 

TTC / 


1/3 


(long, case) (46) 



Fig. 1. Lorentz factor F of the flow as function of radius 
r for the longitudinal and transversal cases. Model pa¬ 
rameters are tro = 100, e = 0.1, p'^ = 0.5, 51 = lO'^s”^, 
1? = 10° and rg = 10^ cm. The solid lines are the numerical 
solutions of (|^) while the dashed lines are the approxi¬ 
mations (^ and (|4^). The vertical lines correspond to 
the photospheric radii: dotted for the transversal case and 
dashed-dotted for the longitudinal case model. The steep 
dashed-dotted line represents the m ~ r law which is ex¬ 
pected in the classic non-magnetic optically thick fireball 
models (Paczyhski 1986). 


and 



(transversal case). (47) 


The function u for the longitudinal case can be de¬ 
scribed as a broken power-law as can be seen in Fig. ^ 
In the domain rg r, mq m <C Mqo, which we have 
regarded anyway for finding the solution, the 4-velocity is 
well approximated by 

3 \ 

— efluoaT ■ (48) 

TTC J 



3.6. The length scale for the acceleration 

In the longitudinal case the dissipation stops approxi¬ 
mately where the rising power-law part of the u functions 
( ^ reaches the Mqo limit ( ^ ) . By using (|3|) and (^ one 
can write down this saturation radius as 

^ ^ TTC gg(l - /r^)^ 

3 ell 3 eQ 

= 7.85-10i2cm-(e_in4)”V3 2 . (49) 


u = 
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The complete analytical approximation for the longitudi¬ 
nal case reads 


u = 


Moo ir/vsr)^^^ for r < 

Moo for r > Tsr 


(50) 


4. Analysis of the results 

4.1. The photospheric radius 


The photosphere is located where the optical depth 
reaches a value of 1. The optical depth depends on the 
density and the the radial velocity u and must be inte¬ 
grated from a finite radius to infinity. Because we only 
want to know the photospheric radius within a factor of, 
say 2, we define it to be where the mean free path of a 
photon equals the distance from the source r. In the co¬ 
moving frame a photon sees the mass density p and the 
mean free path for Thompson scattering is 1/{kp). The 
source distance in this frame is r/T so that the photo¬ 
sphere is located at 


f^Pph 


which yields 
kL 2 

3/2 ^ “ph’^ph ■ 

c3fTo' 


(51) 


(52) 


If we neglect not only the initial velocity but also the 
initial radius ro compared to the photospheric radius, ( ^^ 
and (^) simplify to 

itpij = —ef2(l — p )(Tg' Tph (longitudinal case) , (53) 


7/2 - 

Mph — 


4e(l - 


In 



(transversal case) . 


(54) 


Together with condition (|^) at the photosphere we arrive 
at the equations for the photospheric radius and the 4- 
velocity in the longitudinal case 


'^ph — 


3k 


1 1/5 


= 119- 


I'ph = 


:n(l — p)L 


1- p^\ 

C-iVti ( ^ ^ j T50 

1/5 


]Ji 

9c7 


1/5 


- 3/2 


= 1.05 • 10^^ cm 


e_ 1 H 4 


0.5 


1 -2/5 


^3/5 -3/2 
1^50 ^ 0,2 


(55) 


(56) 


Note that the 4-velocity at the photosphere Wph does not 
depend on the initial Poynting flux ratio ctq and only 
weakly on L. 

For the transversal case the flow velocity always de¬ 
pends greatly on the initial radius rg. The dissipation time 
scale is rtr ~ ru and most energy is released at small r near 


the source. The acceleration depends crucially on the onset 
of the dissipation and therefore on rg. In our simple model 
rg and d are not well determined by physical arguments 
so that the transversal case is rather uncertain and highly 
speculative. One cannot write down robust equations for 
the photosphere like in the longitudinal case without many 
degrees of freedom. 


4.2. Energy available for prompt radiation 

The energy dissipated beyond the photospheric radius is 

Lb = (moo - Mph) Mc^ . (57) 

Using (p^, (|^ and (^ ) this yields 

Ld = e (1 - L (58) 


with 


e = 1 — 


3k eflL 


1/5 


-3/2 


(59) 


For a Poynting flux dominated flow, the magnetic energy 
flux equals the total energy flux L. Of this, a fraction 1—p? 
is dissipated internally while a fraction e(l — p2) can be 
converted to radiation beyond the photosphere. Thus e is 
an efficiency factor, which gives the ratio between the en¬ 
ergy dissipated in the optically thin domain and the total 
dissipated energy. Efficient conversion of free magnetic en¬ 
ergy into non-thermal radiation can happen if e is of the 
order unity which requires that the second term in ( p^ is 
small: 

/ 1 ^ \ 

0.24 • ^ ^ 

or written differently 

/i _ 2\ -8/15 

CTo > 39 • (£- 1114 ^ 50 )^^^^ ( Q ^ j 

where e_i = e/ 0 . 1 , O 4 = n/(l 0 '^s“^), L 50 = 

Lj (l0®°ergs“^ sterad”^) and cro ,2 = uo/lOO are parame¬ 
ters scaled to fiducial GRB values. 

When & is satisfied, part of the magnetic energy is 
released beyond the photosphere, and powers the prompt 
radiation. If it is not satisfied, the energy is released in¬ 
side the photosphere and is converted, instead, into bulk 
kinetic energy. Some other means of conversation into ra¬ 
diation is then needed, such as internal shocks. Since the 
dependence on parameters other than the initial Poynting 
flux ratio ctq is small in (pl|), we conclude that efficient 
powering of prompt radiation by magnetic dissipation in 
GRB is possible for erg ^ 100. 


4.3. Longitudinal and transversal cases in comparison 

The major difference between the longitudinal and the 
transversal case is the different dissipation time scale. 
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Fig. 2. The influence of rg on the Lorentz factor: the solid 
lines correspond to longitudinal case solutions and the 
dashed one to transversal case solutions with 1.5-lO^cm < 
ro < 1.2 • 10® cm. The constants used where erg = 100, 
/i® = 0.5, e = 0.1, n = 10''’s“^. For the transversal cases 
19 = 27rcy^/(Oro) is chosen so that the initial acceler¬ 
ation (slope at ro) is the same as in the corresponding 
longitudinal cases. 


Fig. 3. The Lorentz factor T for different parameters 
a € {0,0.2,0.4,0.6,0.8,1}. The other wind parameters 
were set to cro = 100, y? = 0.5, e = 0.1, ro = 10® cm, 
n = 10^ s“^ (corresponding tod = 10.8° in the transversal 
description). The different winds all start with the same 
dissipation rate so that they show the same initial accel¬ 
eration. Too is marked by a horizontally dotted line. 


While the decay time scale for the longitudinal case (Q) 
is Tio ~ V? and therefore limited by it < itoo, the time 
scale for the transversal case Ttr ~ ru is not limited. At 
small radii it starts at low values but grows then to infin¬ 
ity. This major difference is visualised in Fig. where the 
flow Lorentz factor is plotted depending on the radius. 

As mentioned in Sect. 4.1 the transversal case depends 
strongly on the initial radius rg. This is seen in Fig. ^ 
where the numerical solutions of (^^ are shown for various 
initial radii. While all longitudinal case solutions merge 
toward the u ~ power-law there is a large spread in 
the transversal case solutions. 

The longitudinal and transversal cases are the two lim¬ 
its for a general case where both kinds of scaling of the 
decay time scale occur. One can model the mixing of both 
cases by writing the dissipation time scale as 


r = klLy(}L 

^roj \uo 


2 - 0 . 


(62) 


where 0 < a < 1 is a dimensionless parameter which 
determines the mixing, a = 0 corresponds to the pure 
longitudinal case and a = 1 to the transversal case. The 
constant k can be written depending on the corresponding 
model parameters as 


or as 


k = 


drg,/a^ 

ec 


(64) 


Figure ^ shows the velocity profiles for various a val¬ 
ues. All graphs result from a numerical integration of (^. 
Beyond the photosphere, assuming it is around 10^® cm, 
the dissipation is only efheient if the held variation does 
not point in transversal direction. The efhciency estima¬ 
tion in ( ^^ is therefore an upper limit for the general case 
with a y 0. 


4.4. The validity of the MHD condition 

Because we work with the ideal MHD approximation we 
have to make sure that there are enough charges in the 
flow to make up the required electric current density. 
Because the reconnection processes will destroy the or¬ 
dered initial field configuration quickly it does not make 
much sense to consider this configuration throughout the 
flow. But one can at least estimate needed currents by 
at a sinusoidal wave in the equatorial plane. In 
we derived the limiting radius where the MHD 
condition breaks down by using a constant flow speed and 
assumed fi = 0. The condition that enough charges are 
available to carry the current is 


looking 
Paper 1 


k = 


27070 

efl 


flB 2pec 

dTrF^ nip 


(63) 


(65) 
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This yields the radius up to which the MHD approxima¬ 
tion holds: 


rMHD = 


4e-\/7rL 

mpC^/^O r , 3/2 

P ^ 1 - u/a^' 


.3/2 

(66) 

0 


1^-1 
^0.2 ■ 

(67) 


always fulfilled, as in Paper I 



connection is inefficient for the Crab wind seems to contra¬ 
dict our result, that it efficiently accelerates the GRB out¬ 
flow. The reason for that is the different initial Poynting 
flux values used for the Crab pulsar and in our study, ao is 
the critical parameter controlling the flnal Lorentz factor 
and the spatial size of the accelerating wind. 

Discussing the Crab pulsar wind in detail and specu¬ 
lating why reconnection fails is beyond the scope of the 
present paper. Instead, we simply take the flow parame¬ 
ter values cro,2 = 400, D4 = 0.02 from Lyubarsky & Kirk 
|(200l| ) and show that our model gives basically the same 
result as the stri ped wind model. However, see Lyubarsky 
|fc Eichler (200 1| ) for a critical revision of the Crab pulsar 
wind parameters. Equation yields for the 4-velocity 
at the observed termination shock at r = 3 • 10^^ cm 


-= 0.023 e 


- 1/3 

-1 


— = 920e 

Uo 


- 1/3 

-1 


( 68 ) 


Lyubarsky & Kirk (2001). The different acceleration laws 
of the two models does not change the picture. The ob¬ 
served pulsar wind bubble is to small to allow for efficient 
reconnection. 

The radius Csr from ( ^ ) denotes the radius where the 
Poynting flux conversion ends. Its value scales with the 


third power of crp. Plausible Lorentz factors for GRB winds 
of around 10^-10"^ imply cro,2 ~ 0.2-5 (or larger for /r > 0). 
This lowers by 6 orders of magnitudes compared to 
the Crab wind. Thus Vgr ^ 3 ■ 10^^ cm which is smaller 
than the radius r, « 10^® cm where the flow runs into 


Here, we have used the dependence (22) of the magnetic 
field strength on velocity. In (66) r-MHD is written as a 
function of u and depends implicitly on r. At rp where 
u = uq = it starts at the value (|^) and rises strongly 
until the final velocity u = Uoo = is reached where 
?"MHD diverge to 00. For the GRB parameters assumed 
here, we find tmud ?"sr and the MHD approximation is 


the ambient medium (Piran 1999). The requirement on 
(Tp for the dissipation to take place inside a radius r^, can 
be expressed by ( ]4^ ) which yields an upper limit for the 
initial Poynting flux ratio: 


crp ^ 1100 • (ra,i6e-i^^4)^^^ ^ 


0.5 


-2/3 


(69) 


4.5. Comparison with the striped pulsar wind 

Dissipation of magnetic energy was applie d to the Crab 
pulsar wind by Lyubarsky fc Kirk (2001). Their m odel 
setup included a striped pulsar wind (Coroniti 1990) that 
is the equatorial wind of on inclined rotator. This is quite 
similar to our longitudinal case where all Poynting flux can 
decay so that p = 0. The wind starts with T = y/op and 
reaches Too = ctq^^ as in our model. Due to a difference 
approach to model the reconnection rate they obtain a 


For GRBs there is no size problem as for the Crab wind 
and Poynting flux can be efficiently converted. 

5. Discussion 

We have investigated the effect of dissipation of magnetic 
energy in a GRB outflow on the acceleration of the flow. 
Such dissipation is expected if the flow contains small scale 
changes of direction of the field for example when the flow 
is produced by the the rotation of a non-axisymmetric 
magnetic field. The dissipation is governed by the speed 
of fast reconnection, parameterised in our calculations as 
a fraction e 0.1 of the local Alfven speed in the flow. 

Two possibilities for the field geometry in the outflow 
have been considered: a geometry where the changes in the 
small scale field direction occur along the bulk flow direc¬ 
tion, and a geometry where the field variation is transver¬ 
sal to the flow direction. The first mentioned, longitudinal 
case is expected in the equatorial plane of an inclined ro¬ 


When the wind reaches its termination shock only a small 
fraction of the Poynting flux was converted. Though, due 
to the small amount of mass in the flow large acceleration 
occurs and the Lorentz factor increases by almost 4 or¬ 
ders of magnitude. This is the same result as obtained by 


tator as in the ‘striped’ pulsar wind model of Goronit 
1(19901 ). The second, transversal case can be associated 
with a polar outflow where the field line structure resem¬ 
bles a spiral. In both cases there are MHD instabilities 
(tearing and kink instabilities) which lead to reconnection 
processes. They differ only by the functional form of the 
reconnection time scale. 

We find that in any case the process leads to a strong 
increase of the bulk Lorentz factor of the flow. This ac¬ 
celeration is due to the outward decrease of the magnetic 
pressure resulting from the field decay. At the same time, 
the dissipated energy can be released to large extend in 
the optically thin part of the flow beyond its photosphere, 
and can power most if not all of the prompt emission. This 
provides an alternative to the internal shock model. 

The calculation is done for a stationary wind. Why 
this approximation is valid for highly variable objects like 
GRBs is not obvious. The duration of GRBs t is of the 
order of a few seconds. One can approximate the wind as 
stationary within a source distance ct « 10^^ cm. Thus the 
flow up to the photospheric radius is well described by a 
stationary description. Further out the time dependence 
of a real flow will become more important but that topic 
is beyond the scope of this work. 

The outflow with transversal held variation contains 
some additional complications which does not occur in 
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the longitudinal case. The dissipation time scale is pro¬ 
portional to the source distance. This results in a rapid 
energy dissipation near to the source and the velocity pro¬ 
file depends critically on the radius where the dissipation 
sets in. But this initial radius is hard to estimate from first 
principles. 

We have used the spiral-like field geometry of a po¬ 
lar flow as pictured in Paper I to justify the existence of 
transversal field variations. This field geometry occurs for 
a polar outflow of an axisymmetric rotator. The following 
arguments give reasons why this field geometry is rather 
special and may not be important in a general. The kink 
instability leads to a break-down of the ordered spiral field 
configuration. After some Alfven crossing times the field 
geometry will have changed so that the ‘longitudinal’ dis¬ 
sipation time will become important while the ‘transver¬ 
sal’ time scale grows large and can be neglected. On the 
other hand the rotator may not be perfectly aligned and 
non-axisymmetric field components are also present in the 
polar outflow. So, we probably have always longitudinal 
field variations in the flow so that the findings found in 
our treatment of the ‘longitudinal case’ might be much 
more applicable and general. 

We assume that the thermal energy flux is negligible 
compared to the kinetic and Poynting energy flux. The 
temperature is set to zero which simplifies the treatment 
and allows an analytical integration of the dynamic equa¬ 
tions. Setting the thermal pressure gradient artificially to 
zero might appear to underestimate the acceleration. On 
the other hand the energy equation takes care that all re¬ 
leased magnetic energy shows up in kinetic form. In fact, 
we overestimate the acceleration by doing so because the 
energy part converted into heat reduces the the gain of 
kinetic energy in the flow. Another physical argument ex¬ 
plains why the flow stays cold: The acceleration expressed 
by the scaling of the Lorentz factor gives T ^ for our 
model. The release of magnetic energy must therefore also 
scale with . In contrast to that, purely thermal acceler¬ 
ation by adiabatic cooling leads to more rapid flow acceler¬ 


ation where the Lorentz factor scales like P oc r ( Paczyhski 
1986|)~[rhus, heating proceeds slower than adiabatic cool¬ 
ing so that the thermal pressure gradient is not important 
compared to the magnetic pressure gradient which drives 
the flow. The reason why Lyubarsky fc Kirk (200 1| ) find 
a faster acceleration of T ~ in a similar model lies in 
the different reconnection prescription and is not due to 
their inclusion of thermal pressure. 

In the optically thin regime part of the dissipated en¬ 
ergy radiates away. There, the model over-estimates the 
gain of kinetic energy. We cannot give arguments how 
much dissipated energy escapes as prompt radiation so 
that the total amount of released energy gives only an 
upper limit on the Lorentz factor. 

The photospheric radius determines the lower limit on 
radius for the region in which non-thermal radiation is ex¬ 
pected to originate. For typical GRB parameters describ¬ 
ing the total luminosity, the baryon loading, the fraction 
of dissipatable energy and the reconnection rate one finds 


that a considerable amount of dissipation takes place in 
the optically thin region. Part of the dissipated energy 
is converted into non-thermal radiation. The remainder 
still leads to an acceleration of the flow. This acceleration 
is caused by the magnetic pressure gradient induced by 
the field dissipation. Since the acceleration continues out¬ 
side the photosphere up to the radius where all the free 
magnetic energy is used up this non-thermal radiation is 
emitted from matter with different Lorentz factors. The 
observable spectrum in thus smeared out compared to a 
spectrum from a uniformly moving medium. For a more 
sound analysis of this topic one needs a model for the 
radiation process. 

The Poynting flux conversion happens at radii r ^ 
10 ^® cm which is inside the distance « 10^® cm where the 
GRB outflow is expected to run into the external medium. 
Thus, the Poynting flux can be converted efficiently. But 
by applying the model to the Grab pulsar wind we come 
to the same conclusions as Lyubarsky fc Kirk (2001 ): The 
conversion is inefficient since the observed pulsar wind 
bubble is to small to contain the whole region where re¬ 
connection takes place. For the Grab pulsar the assumed 
initial Poynting flux ratio is larger than for GRBs lead¬ 
ing to a much longer reconnection phase. The presented 
model does not settle this Crab wind problem. 

The most important parameter which controls the 
amount of energy dissipated beyond the photosphere is 
the initial Poynting flux to kinetic energy flux ratio. If its 
value is around 100 or greater much non-thermal, prompt 
emission is produced. If its value is of the order of 10, 
however, all the Poynting flux energy is converted into ki¬ 
netic energy and thermal radiation. Only prompt thermal 
emission and afterglow emission is expected in this case. 
The initial Poynting flux ratio is a measure for the baryon 
loading in a sense that a high baryon loading corresponds 
to a low initial Poynting flux ratio. Observations indicate 
that X-ray flashes and X-ray rich GRBs are very simi¬ 
lar phenomena whi ch probably differ only by the amount 
of baryon loading ( Beise et al. 2001 ). In the context of 
our model, X-ray flashes can be associated with low ini¬ 
tial Poynting flux ratios. In this case, the X-ray emission 
is thermal radiation from the photosphere. Increasing the 
initial Poynting flux ratio leads to the emission of non- 
thermal 7-rays in the optically thin region, thus produc¬ 
ing X-ray rich and regular GRBs. If afterglows of X-ray 
flashes could be observed they would yield information 
about the connection to regular GRBs. Afterglows depend 
less strongly on the initial Poynting flux ratio but rather 
on the total luminosity of the outflow. Thus, X-ray flash 
afterglows should be similar to afterglows of regular GRBs 
according to our model. In a future work we will investi¬ 
gate the thermal emission more quantitatively. 

The model predicts black-body radiation originating 
from the photosphere of the flow. We can calculate the 
radius of the photosphere and the Lorentz factor of the 
flow there. Together with the temperature one is able to 
calculate the luminosity if the thermal radiation. Since our 
approximation treats the flow as cold we cannot give quan- 
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titative results in this respect. Though, one finds that the 
Lorentz factor at the photosphere depends only weakly on 
the model parameters. Therefore, the observable tempera¬ 
ture fcTobs = /{1 + z) of the thermal component of a 

GRB depends primarily on the redshift z and the temper¬ 
ature in the comoving frame T. This result simplifies the 
task to disentangle the effects of different model parame¬ 
ters on the temperature. A detailed, quantitative analysis 
of the thermal radiation will be done in a following study. 

Acknowledgements. I thank H.G. Spruit for enlightening dis¬ 
cussions and the critical reading of the manuscript. 
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